Nonlinear Enzyme Kinetics Can Lead to High Metabolic Flux Control Coefficients: Implications for the Evolution of Dominance

In a classic study, Kacser & Burns (1981, Genetics 97, 639-666) demonstrated that given certain plausible assumptions, the flux in a metabolic pathway was more or less indifferent to the activity of any of the enzymes in the pathway taken singly. It was inferred from this that the observed dominance of most wild-type alleles with respect to loss-of-function mutations did not require an adaptive, meaning selectionist, explanation. Cornish-Bowden (1987, J. theor. Biol. 125, 333-338) showed that the Kacser-Burns inference was not valid when substrate concentrations were large relative to the relevant Michaelis constants. We find that in a randomly constructed functional pathway, even when substrate levels are small, one can expect high values of control coefficients for metabolic flux in the presence of significant nonlinearities as exemplified by enzymes with Hill coefficients ranging from two to six, or by the existence of oscillatory loops. Under these conditions the flux can be quite sensitive to changes in enzyme activity as might be caused by inactivating one of the two alleles in a diploid. Therefore, the phenomenon of dominance cannot be a trivial ''default'' consequence of physiology but must be intimately linked to the manner in which metabolic networks have been moulded by natural selection.

trans gene regulation in adaptive evolution: A genetic algorithm model

This is a continuation of earlier studies on the evolution of infinite populations of haploid genotypes within a genetic algorithm framework. We had previously explored the evolutionary consequences of the existence of indeterminate—“plastic”—loci, where a plastic locus had a finite probability in each generation of functioning (being switched “on”) or not functioning (being switched “off”). The relative probabilities of the two outcomes were assigned on a stochastic basis. The present paper examines what happens when the transition probabilities are biased by the presence of regulatory genes. We find that under certain conditions regulatory genes can improve the adaptation of the population and speed up the rate of evolution (on occasion at the cost of lowering the degree of adaptation). Also, the existence of regulatory loci potentiates selection in favour of plasticity. There is a synergistic effect of regulatory genes on plastic alleles: the frequency of such alleles increases when regulatory loci are present. Thus, phenotypic selection alone can be a potentiating factor in a favour of better adaptation.

Genetic studies of the PRP17 gene of Saccharomyces cerevisiae: a domain essential for function maps to a nonconserved region of the protein

The PRP17 gene product is required for the second step of pre-mRNA splicing reactions. The C-terminal half of this protein bears four repeat units with homology to the beta transducin repeat. Missense mutations in three temperature-sensitive prp17 mutants map to a region in the N-terminal half of the protein. We have generated, in vitro, 11 missense alleles at the beta transducin repeat units and find that only one affects function in vivo. A phenotypically silent missense allele at the fourth repeat unit enhances the slow-growing phenotype conferred by an allele at the third repeat, suggesting an interaction between these domains. Although many missense mutations in highly conserved amino acids lack phenotypic effects, deletion analysis suggests an essential role for these units. Only mutations in the N-terminal nonconserved domain of PRP17 are synthetically lethal in combination with mutations in PRP16 and PRP18, two other gene products required for the second splicing reaction. A mutually allele-specific interaction between Prp17 and snr7, with mutations in U5 snRNA, was observed. We therefore suggest that the functional region of Prp17p that interacts with Prp18p, Prp16p, and U5 snRNA is the N terminal region of the protein.

Nucleotide sequence and expression inE. coli of the complete P4 type VP4 from a G2 serotype human rotavirus

The complete sequence of a P4 type VP4 gene from a G2 serotype human rotavirus, IS2, isolated in India has been determined. Although the IS2 VP4 is highly homologous to the other P4 type alleles, it contained acidic amino acid substitutions at several positions that make it acidic among the P4 type alleles that are basic. Moreover, comparative sequence analysis revealed unusual polymorphism in members of the P4 type at amino acid position 393 which is highly conserved in members of other VP4 types. To date, expression of complete VP4 inE. coli has not been achieved. In this study we present successful expression inE. coli of the complete VP4 as well as VP8* and VP5* cleavage subunits in soluble form as fusion proteins of the maltose-binding protein (MBP) and their purification by single-step affinity chromatography. The hemagglutinating activity exhibited by the recombinant protein was specifically inhibited by the antiserum raised against it. Availability of pure VP4 proteins should facilitate development of polyclonal and monoclonal antibodies (MAbs) for P serotyping of rotaviruses.

Mitochondrial and nuclear markers suggest Hanuman langur (Primates: Colobinae) polyphyly: Implications for their species status

Recent molecular studies on langurs of the Indian subcontinent suggest that the widely-distributed and morphologically variable Hanuman langurs (Semnopithecus entellus) are polyphyletic with respect to Nilgiri and urple-faced langurs. To further investigate this scenario, we have analyzed additional sequences of mitochondrial cytochrome b as well as nuclear protamine P1 genes from these species. The results confirm Hanuman langur polyphyly in the mitochondrial tree and the nuclear markers suggest that the Hanuman langurs share protamine P1 alleles with Nilgiri and purple-faced langurs. We recommend provisional splitting of the so-called Hanuman langurs into three species such that the taxonomy is consistent with their evolutionary relationships.

Analysis of the impact of a uracil DNA glycosylase attenuated in AP-DNA binding in maintenance of the genomic integrity in Escherichia coli

Uracil DNA glycosylase (Ung)initiates the uracil excision repair pathway. We have earlier characterized the Y66W and Y66H mutants of Ung and shown that they are compromised by similar to 7- and similar to 170-fold, respectively in their uracil excision activities. In this study, fluorescence anisotropy measurements show that compared with the wild-type, the Y66W protein is moderately compromised and attenuated in binding to AP-DNA. Allelic exchange of ung in Escherichia coli with ung::kan, ungY66H:amp or ungY66W:amp alleles showed similar to 5-, similar to 3.0- and similar to 2.0-fold, respectively increase in mutation frequencies. Analysis of mutations in the rifampicin resistance determining region of rpoB revealed that the Y66W allele resulted in an increase in A to G (or T to C) mutations. However, the increase in A to G mutations was mitigated upon expression of wild-type Ung from a plasmid borne gene. Biochemical and computational analyses showed that the Y66W mutant maintains strict specificity for uracil excision from DNA. Interestingly, a strain deficient in AP-endonucleases also showed an increase in A to G mutations. We discuss these findings in the context of a proposal that the residency of DNA glycosylase(s) onto the AP-sites they generate shields them until recruitment of AP-endonucleases for further repair.

Confirmation of linkage and refinement of the RP28 locus for autosomal recessive retinitis pigmentosa on chromosome 2p14-p15 in an Indian family

PURPOSE: To report the linkage analysis of retinitis pigmentosa (RP) in an Indian family. METHODS: Individuals were examined for symptoms of retinitis pigmentosa and their blood samples were withdrawn for genetic analysis. The disorder was tested for linkage to known 14 adRP and 22 arRP loci using microsatellite markers. RESULTS: Seventeen individuals including seven affecteds participated in the study. All affected individuals had typical RP. The age of onset of the disease ranged from 8-18 years. The disorder in this family segregated either as an autosomal recessive trait with pseudodominance or an autosomal dominant trait. Linkage to an autosomal recessive locus RP28 on chromosome 2p14-p15 was positive with a maximum two-point lod score of 3.96 at theta=0 for D2S380. All affected individuals were homozygous for alleles at D2S2320, D2S2397, D2S380, and D2S136. Recombination events placed the minimum critical region (MCR) for the RP28 gene in a 1.06 cM region between D2S2225 and D2S296. CONCLUSIONS : The present data confirmed linkage of arRP to the RP28 locus in a second Indian family. The RP28 locus was previously mapped to a 16 cM region between D2S1337 and D2S286 in a single Indian family. Haplotype analysis in this family has further narrowed the MCR for the RP28 locus to a 1.06 cM region between D2S2225 and D2S296. Of 15 genes reported in the MCR, 14 genes (KIAA0903, OTX1, MDH1, UGP2, VPS54, PELI1, HSPC159, FLJ20080, TRIP-Br2, SLC1A4, KIAA0582, RAB1A, ACTR2, and SPRED2) are either expressed in the eye or retina. Further study needs to be done to test which of these genes is mutated in patients with RP linked to the RP28 locus.

An automated framework for understanding structural variations in the binding grooves of MHC class II molecules

Background: MHC/HLA class II molecules are important components of the immune system and play a critical role in processes such as phagocytosis. Understanding peptide recognition properties of the hundreds of MHC class II alleles is essential to appreciate determinants of antigenicity and ultimately to predict epitopes. While there are several methods for epitope prediction, each differing in their success rates, there are no reports so far in the literature to systematically characterize the binding sites at the structural level and infer recognition profiles from them. Results: Here we report a new approach to compare the binding sites of MHC class II molecules using their three dimensional structures. We use a specifically tuned version of our recent algorithm, PocketMatch. We show that our methodology is useful for classification of MHC class II molecules based on similarities or differences among their binding sites. A new module has been used to define binding sites in MHC molecules. Comparison of binding sites of 103 MHC molecules, both at the whole groove and individual sub-pocket levels has been carried out, and their clustering patterns analyzed. While clusters largely agree with serotypic classification, deviations from it and several new insights are obtained from our study. We also present how differences in sub-pockets of molecules associated with a pair of autoimmune diseases, narcolepsy and rheumatoid arthritis, were captured by PocketMatch(13). Conclusion: The systematic framework for understanding structuralvariations in MHC class II molecules enables large scale comparison of binding grooves and sub-pockets, which is likely to have direct implications towards predicting epitopes and understanding peptide binding preferences.

An automated framework for understanding structural variations in the binding grooves of MHC class II molecules

Background: MHC/HLA class II molecules are important components of the immune system and play a critical role in processes such as phagocytosis. Understanding peptide recognition properties of the hundreds of MHC class II alleles is essential to appreciate determinants of antigenicity and ultimately to predict epitopes. While there are several methods for epitope prediction, each differing in their success rates, there are no reports so far in the literature to systematically characterize the binding sites at the structural level and infer recognition profiles from them. Results: Here we report a new approach to compare the binding sites of MHC class II molecules using their three dimensional structures. We use a specifically tuned version of our recent algorithm, PocketMatch. We show that our methodology is useful for classification of MHC class II molecules based on similarities or differences among their binding sites. A new module has been used to define binding sites in MHC molecules. Comparison of binding sites of 103 MHC molecules, both at the whole groove and individual sub-pocket levels has been carried out, and their clustering patterns analyzed. While clusters largely agree with serotypic classification, deviations from it and several new insights are obtained from our study. We also present how differences in sub-pockets of molecules associated with a pair of autoimmune diseases, narcolepsy and rheumatoid arthritis, were captured by PocketMatch(13). Conclusion: The systematic framework for understanding structural variations in MHC class II molecules enables large scale comparison of binding grooves and sub-pockets, which is likely to have direct implications towards predicting epitopes and understanding peptide binding preferences.

Molecular analysis of Friedreich's ataxia locus in the Indian population

Friedreich's ataxia (FRDA) is an autosomal recessive neurodegenerative disorder caused by expansion of GAA repeats in the frataxin gene. We have carried out the first molecular analysis at the Friedreich's ataxia locus in the Indian population. Materials and methods - Three families clinically diagnosed for Friedreich's ataxia were analyzed for GAA expansion at the FRDA locus. The distribution of GAA repeats was also estimated in normal individuals of Indian origin. Results - All patients clinically diagnosed for Friedreich's ataxia were found to be homozygous for GAA repeat expansion. The GAA repeat in the normal population show a bimodal distribution with 94% of alleles ranging from 7-16 repeats. Conclusion - Indian patients with expansion at the FRDA locus showed typical clinical features of Friedreich's ataxia. The low frequency of large normal alleles (6%) could indicate that the prevalence of this disease in the Indian population is likely to be low.

The evolution of genomic imprinting

We explore three possible pathways for the evolution of genomic imprinting. (1) Imprinting may be advantageous in itself when imprinted and unimprinted alleles of a locus confer different phenotypes. If a segment of DNA is imprinted in the gametes of one sex but not in those of the other, it might lead to effects correlated with sexual dimorphism. More fundamentally, in certain organisms, sex determination might have evolved because of imprinting. When imprinting leads to chromosome elimination or inactivation and occurs in some embryos but not in others, two classes of embryos, differing in the number of functional gene copies, would result. A model for sex determination based on inequality in the actual or effective copy-number of particular noncoding, regulatory sequences of DNA has been proposed (Chandra, Proc. natn. Acad. Sci. U.S.A. 82. 1165–1169 and 6947–6949, 1985). Maternal control of offspring sex is another possible consequence of imprinting; this would indicate a potential role for imprinting in sex ratio evolution. (2) Genes responsible for imprinting may have pleiotropic effects and they may have been selected for reasons other than their imprinting ability. Lack of evidence precludes further consideration of this possibility. (3) Imprinting could have co-evolved with other traits. For instance, gamete-specific imprinting could lead to a lowered fitness of androgenetic or gynogenetic diploids relative to the fitness of ‘normal’ diploids. This in turn would reinforce the evolution of anisogamy. The reversibility of imprinting raises the possibility of occasional incomplete or improper erasure. If the site of imprinting is the egg – as appears to be the case with the human X (Chandra and Brown, Nature 253. 165–168, 1975) – either improper imprinting or improper erasure could lead to unusual patterns of inheritance (as in the fragile-X syndrome) or fitness effects skipping generations.

Comparative Nucleotide and Amino Acid Sequence Analysis of the Sequence-Specific RNA-Binding Rotavirus Nonstructural Protein NSP3

NSP3, an acidic nonstructural protein, encoded by gene 7 has been implicated as the key player in the assembly of the 11 viral plus-strand RNAs into the early replication intermediates during rotavirus morphogenesis. To date, the sequence or NSP3 from only three animal rotaviruses (SA11, SA114F, and bovine UK) has been determined and that from a human strain has not been reported. To determine the genetic diversity among gene 7 alleles from group A rotaviruses, the nucleotide sequence of the NSP3 gene from 13 strains belonging to nine different G serotypes, from both humans and animals, has been determined. Based on the amino acid sequence identity as well as phylogenetic analysis, NSP3 from group A rotaviruses falls into three evolutionarily related groups, i.e., the SA11 group, the Wa group, and the S2 group. The SA 11/SA114F gene appears to have a distant ancestral origin from that of the others and codes for a polypeptide of 315 amino acids (aa) in length. NSP3 from all other group A rotaviruses is only 313 aa in length because of a 2-amino-acid deletion near the carboxy-terminus, While the SA114F gene has the longest 3' untranslated region (UTR) of 132 nucleotides, that from other strains suffered deletions of varying lengths at two positions downstream of the translational termination codon. In spite of the divergence of the nucleotide (nt) sequence in the protein coding region, a stretch of about 80 nt in the 3' UTR is highly conserved in the NSP3 gene from all the strains. This conserved sequence in the 3' UTR might play an important role in the regulation of expression of the NSP3 gene. (C) 1995 Academic Press, Inc.

Sibling rivalry between seeds within a fruit: Some population genetic models

Competition between seeds within a fruit for parental resources is described using one-locus-two-allele models. While a �normal� allele leads to an equitable distribution of resources between seeds (a situation which also corresponds to the parental optimum), the �selfish� allele is assumed to cause the seed carrying it to usurp a higher proportion of the resources. The outcome of competition between �selfish� alleles is also assumed to lead to an asymmetric distribution of resources, the �winner� being chosen randomly. Conditions for the spread of an initially rare selfish allele and the optimal resource allocation corresponding to the evolutionarily stable strategy, derived for species with n-seeded fruits, are in accordance with expectations based on Hamilton�s inclusive fitness criteria. Competition between seeds is seen to be most intense when there are only two seeds, and decreases with increasing number of seeds, suggesting that two-seeded fruits would be rarer than one-seeded or many-seeded ones. Available data from a large number of plant species are consistent with this prediction of the model.

The carboxy terminal WD domain of the pre-mRNA splicing factor Prp17p is critical for function

In Saccharomyces cerevisiae, Prp17p is required for the efficient completion of the second step of pre-mRNA splicing. The function and interacting factors for this protein have not been elucidated. We have performed a mutational analysis of yPrp17p to identify protein domains critical for function. A series of deletions were made throughout the region spanning the N-terminal 158 amino acids of the protein, which do not contain any identified structural motifs. The C-terminal portion (amino acids 160–455) contains a WD domain containing seven WD repeats. We determined that a minimal functional Prp17p consists of the WD domain and 40 amino acids N-terminal to it. We generated a three-dimensional model of the WD repeats in Prp17p based on the crystal structure of the [beta]-transducin WD domain. This model was used to identify potentially important amino acids for in vivo functional characterization. Through analysis of mutations in four different loops of Prp17p that lie between [beta] strands in the WD repeats, we have identified four amino acids, 235TETG238, that are critical for function. These amino acids are predicted to be surface exposed and may be involved in interactions that are important for splicing. Temperature-sensitive prp17 alleles with mutations of these four amino acids are defective for the second step of splicing and are synthetically lethal with a U5 snRNA loop I mutation, which is also required for the second step of splicing. These data reinforce the functional significance of this region within the WD domain of Prp17p in the second step of splicing.

An analysis of substitution, deletion and insertion mutations in cancer genes

Cancer-associated mutations in cancer genes constitute a diverse set of mutations associated with the disease. To gain insight into features of the set, substitution, deletion and insertion mutations were analysed at the nucleotide level, from the COSMIC database. The most frequent substitutions were c -> t, g -> a, g -> t, and the most frequent codon changes were to termination codons. Deletions more than insertions, FS (frameshift) indels more than I-F (in-frame) ones, and single-nucleotide indels, were frequent. FS indels cause loss of significant fractions of proteins. The 5'-cut in FS deletions, and 5'-ligation in FS insertions, often occur between pairs of identical bases. Interestingly, the cut-site and 3'-ligation in insertions, and 3'-cut and join-pair in deletions, were each found to be the same significantly often (p < 0.001). It is suggested that these features aid the incorporation of indel mutations. Tumor suppressors undergo larger numbers of mutations, especially disruptive ones, over the entire protein length, to inactivate two alleles. Proto-oncogenes undergo fewer, less-disruptive mutations, in selected protein regions, to activate a single allele. Finally, catalogues, in ranked order, of genes mutated in each cancer, and cancers in which each gene is mutated, were created. The study highlights the nucleotide level preferences and disruptive nature of cancer mutations.